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Context. In a flat space, it has been shown heuristically that the global topology of comoving space can affect the 
dynamics expected in the weak-field Newtonian limit, inducing a weak acceleration effect similar to dark energy. 
Aims. Does a similar effect occur in the case of the Poincare dodecahedral space, which is a candidate model of comoving 
space for solving the missing fluctuations problem observed in cosmic microwave background all-sky maps? Moreover, 
does the effect distinguish the Poincare space from other well-proportioned spaces? 

Methods. The acceleration effect in the Poincare space S 3 /I* is studied, using a massive particle and a nearby test 
particle of negligible mass. Calculations are made in S 3 embedded in R 4 . The weak-limit gravitational attraction on 
a test particle at distance r is set to be oc [Rc s'm(r / Rc)]~ 2 rather than oc r -2 , where Rc is the curvature radius, in 
order to satisfy Stokes' theorem. A finite particle horizon large enough to include the adjacent topological images of 
the massive particle is assumed. The regular, flat, 3-torus T s is re-examined, and two other well-proportioned spaces, 
the octahedral space S 3 /T* , and the truncated cube space s 3 /o*, are also studied. 

Results. The residual gravity effect occurs in all four cases. In a perfectly regular 3-torus of side length L a , and in 
the octahedral and truncated cube spaces, the highest order term in the residual acceleration is the third-order term 
in the Taylor expansion in powers of r/L a (3-torus), or r/Rc, respectively. However, the Poincare dodecahedral space 
is unique among the four spaces. The third order cancels, leaving the fifth order term ~ ±300(r/i?c) 5 as the most 
significant. 

Conclusions. Not only are three of the four perfectly regular well-proportioned spaces better balanced than most other 
multiply connected spaces in terms of the residual gravity acceleration effect by a factor of about a million (setting 
r/L a = r/Rc ~ 10~ 3 ), but the fourth of these spaces is about ten thousand times better balanced than the other 
three. This is the Poincare dodecahedral space. Is this unique dynamical property of the Poincare space a clue towards 
a theory of cosmic topology? 
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ual gravity effect might contribute an additional criterion 
for comparing different manifolds. 

Here, the residual gravity effect in the octahedral space, 
the truncated cube space, and the Poincare space is studied 
by considering the dynamics of a negligible-mass test par- 
ticle near a massive particle. As in Roukema et al. (2007), 
this approach can be considered as an heuristic model for 
a positively curved space that is homogeneous except for 
a small neighbourhood around one point, in which a posi- 
tive density fluctuation of matter has collapsed into a high- 
density, nearly point-like object, in excess of the underlying 
homogeneous density distribution. To satisfy Stokes' theo- 
rem, the weak-limit gravitational attraction in the three 
spherical spaces is set to be proportional to 



(i) 



rather than r~ 2 (see Sect. I2.1[) . where Rc is the curva- 
ture radius. Numerical calculations are made in the uni- 
versal covering space S 3 embedded in R for convenience. 
Accelerations are calculated in the tangent 3-space at the 
position of the test particle, represented as 4- vectors in R 4 . 
The T 3 model discussed in Roukema et al. (2007) is revis- 
ited, using the covering space R 3 . 

The method is described in more detail in Sect. [U 
Residual accelerations estimated to the third order in the 
Taylor expansion of the fractional displacement are pre- 
sented for the 3-torus in Sect. 13. f I Residual accelerations 
to fifth order for individual pairs of opposite images in 
the spherical spaces are presented analytically in Sect. 13.21 
Residual accelerations to fifth order for the full set of adja- 
cent images in the spherical spaces are presented numer- 
ically in Sect. 13.31 For the Poincare space, an analytical 
derivation of the residual acceleration to fifth order is also 
given in Sect. 13.3.11 A statistical summary is presented in 
Sect. 13X21 

Discussion and conclusions are given in Sects. 2] and 
[5] respectively. Discussion of spherical, multiply connected 
spaces is available in Weeks (2001), Gausmann et al. (2001), 
Lehoucq et al. (2002), and Riazuelo et al. (2004), and we 
refer the reader to references therein for introductions to 
cosmic topology. Distances are calculated in a spatial sec- 
tion at constant cosmological time in the universal covering 
space (S 3 for the spherical spaces, R 3 for T 3 ), by default 
as "physical" distances, i.e., r = ax, where a(t) is the scale 
factor and x is a comoving distance (i.e. \ is a "proper 
distance" at the present epoch, (Weinberg 1972), equiva- 
lent to "conformal time" if c — 1). Distance units are nor- 
mally presented here in h~ l Mpc or hr 1 Gpc, where the 
Hubble constant is written Hq = 100/ikms _1 /Mpc. The 
fundamental domains of all the spaces are assumed to be 
perfectly regular. The Newtonian gravitational constant is 
written as G. 



2. Method 

Although dynamics in the Poincare dodecahedral space 
S 3 /I* could, in principle, be studied by applying boundary 
conditions that identify opposite faces of a dodecahedral 



fundamental domain of a positive curvature radius Rc with 
one another, this would probably be extremely complicated 
for both analytical and numerical calculations. It is simpler 
to work in the universal covering space S 3 represented as a 
subspace of R 4 . 

As described in Sect. 2.1 of Roukema et al. (2007), sev- 
eral assumptions are required for this heuristic approach. 
Here, the corresponding assumptions are as follows: 

(1) the flat Newtonian approximation of gravity is replaced 
by the equivalent in positively curved space, as de- 
scribed in Eq. (fT]) above; 

(2) , (3) the covering space is S 3 , which is not flat; 

(4) a finite particle horizon i?hor that is just large enough 
to include all adjacent topological images is assumed; 

(5) , (6) identical to those in Roukema et al. (2007): the 
metric is assumed to be that for a perfectly homoge- 
neous model of the same curvature, except that the dis- 
tant, multiple copies of the local massive particle are 
considered to provide a contribution to the local gravi- 
tational potential that may not fully cancel; 

(7), (8) are not needed here since the one-body problem in 
S 3 (Sect. 3.1.1, Roukema et al.) is divergent and only 
briefly mentioned in Sect. 12.11 

2.1. Weak-limit gravity and divergences 

As stated in Eq.fl]), to satisfy Stokes' theorem, the weak- 
limit gravitational attraction towards a single massive par- 
ticle is set to be proportional to [Rc sin(r/i?c)]~ 2 rather 
than r~ 2 . This can be understood as follows. 

Let S 3 be represented by a spherical coordinate system 
centred on the massive particle, so that the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) line element is writ- 
ten 

ds 2 = ~c 2 dt 2 +dr 2 + [i? c sin(r/i? c )] 2 (d6> 2 + cos 2 6>d0 2 ), 

(2) 

where the scale factor a(t) is subsumed into r and Rc, 
since it is not needed explicitly. We can write the weak- 
limit gravitational acceleration towards the massive particle 
as the vector field 



/Mr 



(3) 



by symmetry, where f is the unit radial vector in the tan- 
gent space at a radial "physical" distance r from the origin 
along a spatial geodesic, and 



-f(r) r 2 

lim = 1 (4) 

r— >0+ Gm 

in order to agree with the flat-space Newtonian limit to- 
wards the massive particle of mass m. Let V be the interior 
of a 2-sphere of radius r centred at the massive particle. 
Applying Stokes' theorem to f(r) f on V, which is a sub- 
manifold of S 3 with boundary dV, gives 

div[/(r') r]dV = / /(/) f .dS 

JdV 

f(r')[R c sm(r'/Rc)} 2 d6 cosfld^, 

dV 

f(r)[R c sm{r/Rc)} 2 [ 60 cos 6d(j>, 

JdV 

4TT[R c sm(r/Rc)} 2 f(r), (5) 



1 "Adjacent" is used here to refer to images in copies of the 
fundamental domain that share a face with the "original" copy 
of the fundamental domain. 
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using the FLRW metric. Since the only source of the vector 
field f(r) f is the single massive particle, the left-hand side 
of this equation equals its limit for arbitrarily small (non- 
zero) volumes, i.e. 

/ div[/(r') f] dV = lim / div[/(r') f] dV 

Jy r^0+ Jy 



= lim 47r[i? c sin(Y/i?c)] 2 f(r) 
= -47rGm, 



(0) 



using Eqs (O and (fj| respectively. Equating the right-hand 
sides of Eqs ([5]) and ([6j) gives the weak limit acceleration 
induced by a single massive particle in a positively curved 
space of radius Rc 



r = f(r) r = — 



Gm f 



[Rc sin(r/i? c )] 2 



(7) 



For instantaneous transmission of gravitational signals in 
an infinitely old, static, simply connected space, this implies 
two divergence problems. 

Firstly, a single massive point particle of mass m in a 
perfectly uniform and otherwise empty S 3 yields an infinite 
repulsive force at the antipode of the massive particle, since 
sin(7r) = 0. More generally, an infinite repulsive (attractive) 
force exists for gravitational signals that have travelled odd 
(even) values of j times the half-circumference ttRc from 
the massive object, i.e. j/2 times around the whole space, 
where j S 1. 

Secondly, even if we ignore the model of a zero-size point 
particle as an idealised fiction, a problem remains for a 
negligible-mass test particle near a massive particle. The 
test particle also experiences accelerations from signals that 
have travelled j times around S in the two different direc- 
tions along the great circle passing through the test particle 
and the massive particle, where j e Z. This second diver- 
gence problem is similar to the divergence problem in flat, 
simply connected, infinitely sized, infinitely old, Newtonian 
space with instantaneous transmission of the gravitational 
signal. The total acceleration towards the massive particle 
would be 



's 3 



sin 



r + 2jirRc 



Rc 



- Gm [ R c 

sm 



2j7ri? c \ ' 



Rc J 



Gm 




E 

j=— oo 



(8) 



which is clearly divergent. 

The former divergence, i.e., the j~xRc singularity related 
to the zero size of the massive point particle and the nature 
of positive curvature, is clearly unphysical if we are inter- 
ested in a smoothing length on the scale of a galaxy clus- 
ter. The latter divergence, i.e. that in Eq. ([5]), is removed 
by assumption (4) above. Assumption (4) also removes the 
j7ri?c divergence. A finite particle horizon i?hor just large 
enough to include the adjacent topological images is neces- 
sary for a residual gravity effect to occur, but it does not 
need to be as large as 7ri?c- The adjacent topological im- 
ages in S 3 /T*, S 3 /Q*, and S 3 /I* are at (tt/3)R C , (tt/ 4 )#c 



and (7r/3)i?C: and (i"/5)i?c respectively. Moreover, obser- 
vationally realistic estimates of the total density parame- 
teiQ are consistent with the range that is empirically in- 
teresting for the Poincare space, 1.02 >, f7 to t ^ 1.01. The 
range 1.02 > r2 to t~ 1-01 corresponds to a curvature radius 
of 20h~ 1 Gpc^i?c~30/i _1 Gpc respectively. Hence, for 
the three well-proportioned, spherical spaces, the appro- 
priate horizon distances for the arrival of gravitational sig- 
nals from the adjacent topological images are in the range 
10/i -1 Gpc < i?hor ~ 30/i _1 Gpc, i.e. up to a few times the 
distance to the surface of last scattering tsls ~ 10/i -1 Gpc. 
As noted in Roukema et al. (2007), a moderate amount of 
inflation in the early Universe could be one way of satisfying 
assumption (4). 

2.2. Residual gravity 

The weak-limit residual acceleration on a test particle at r 
near the massive particle of mass m at can be written 



E 

j=l,d j= ax[r,gj{Q)] 



Gi 



Rc sin -j^- 

Rc 



L 3> 



(9) 



where gj — 1, . . . , N are the N = 8, 14, or 12 holonomies 
that map an arbitrary object in S 3 /T*, S 3 fO*, or S 3 /I* re- 
spectively to its adjacent topological images, a is the scale 
factor, x( a i t>) is th e comoving distance between two posi- 
tions a, b in the covering space S 3 , Tj is the unit vector in 
the tangent space to S 3 at r pointing along the geodesic 
towards <?j(0), and f is the acceleration vector in the same 
tangent space. The appendices in Gausmann et al. (2001) 
provide a convenient procedure for calculating this, where 
a 3-manifold S 3 /T and holonomy group T are modelled as 
a tiling and isometry group of M 4 using quaternions and 4- 
dimensional matrices. The tangent space at r is a 3-plane, 
modelled using vectors in M 4 . The unit tangent vector to- 
wards the j-th topological image of the massive particle can 
be written 



f = &-(())- ( ffj - (0).r)r 
J ||a(0)-(S0).r)r|| 

where x = x/||x||. A useful identity is 

-2 

9M-t]})' 



(10) 



(sin-j^- ] = fsinl 
V Rc) VI 



acos 



1 



(11) 



The curvature is set to Rc = 25/i _1 Gpc below unless 
stated otherwise. The horizon i?hor is set sufficiently large 
to include the full set of adjacent images of the massive 
particle, i.e., values of i?hor = l.l(7r/3)i?c, l.l(7r/3)i?c, 
and 1.1(^/5)^0 are used in Eq. © for S 3 /T*, S 3 fO*, and 
S 3 /T, respectively. 



2 For example, f2 to t — 1.010 
Hubble Space Telescope Hq key project data, fitot 



° ° ll, from WMAP 3-year and 
1.015+ 020 



0.016' 

from WMAP 3-year data and Supernova Legacy Survey super- 
novae Type la data (Spergel et al. 2007); fitot = LOSt^ 4 , 
from WMAP 3-year and Arcminute Cosmology Bolometer Array 
Receiver ACBAR data (Reichardt et al. 2008). 
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2.3. Numerical precision 

For the full set of adjacent topological images in the spher- 
ical spaces, Eq. @ is evaluated primarily by a numerical 
approach, using double precision floating-point operations 
where sufficient, and arbitrary precision floating-point op- 
erations where necessary. Typical scales of interest, i.e., for 
typical test particles in a void of large-scale structure, are a 
few tens of megaparsecs. This gives 10~ 4 < r/Rc < 1 CP 2 as 
a useful domain for finding the highest terms in the Taylor 
expansion of the residual gravity for r <§; Rq, where i?c = 
25h~ 1 Gpc as stated above. By setting Rq = G = m = 1 
internally in numerical calculations, nearly equal and op- 
posite accelerations from nearly opposed topological im- 
ages are each of approximately unity order, the fourth and 
fifth order terms at the lower limit r/Rc = fCP 4 should 
therefore be above the numerical noise limit, if the calcu- 
lating precision is well below ~ 10 -16 and 10 -20 , respec- 
tively. These would require precision in the significance well 
above log 2 10 16 « 53 and log 2 f 20 ss 66 bits respectively. 
According to the IEEE 754-1 985 standard, double-precision 
floating-point numbers have 53-bit precision in the signifi- 
cand (including one implicit bit). This would at best pro- 
vide just one bit of information for a fourth order term and 
no information for a fifth order term for r/Rc = 10~ 4 . For 
this reason, an arbitrary precision library is used here to 
examine higher order terms. 

Analytical calculations are also made for the Poincare 
space. We also revisit the T 3 calculation made in Sect. 3.2 
of Roukema et al. (2007), in order to consider higher or- 
der terms. The first order term was found to cancel for a 
perfectly regular T 3 model. 



3. Results 

3.1. T 3 revisited 

In Sect. 3.2 of Roukema et al. (2007), an analytical esti- 
mate of the acceleration from all six adjacent images in a 
T 3 model was shown to be zero to first order when the 
three fundamental lengths of the model are exactly equal. 
However, the lower scatterings of points in Figs. 7 and 8 of 
that paper do not represent numerical error. Recalculation 
of Eqs (11) to (15) of Roukema et al. (2007) to higher 
order shows that for perfectly equal fundamental lengths 
L a = L c = L u , higher order terms in Eq. (15) do not all 
cancel. The lowest order non- cancelling terms are the third 
order terms 

(" : topo)(i,j,fc)e{(±i,o,o),(o,±i,o),(o,o,±i)} 
7Gm 

— £2 •'■ 



[2e 3 x -3e x (e 2 +e 2 z ),2e 3 y -3e y (e 2 x + e 2 z ), 



2^-3e,(e 2 +e 2 )], 



(12) 



where the massive particle is at the origin (0,0,0) in the 
covering space M. 3 and the test particle is slightly displaced 
at {x — e x L ai y — e y L ai z — e z L a ). For test particles dis- 
tributed isotropically around the massive particle, numer- 
ical integration of these third order terms gives the distri- 



3 The IEEE 754-2008 standard recommends the term "signif- 
icand" rather than "mantissa". 
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Fig. 1. Residual acceleration (radial component) induced on a 
test particle near a massive particle such as a cluster of galaxies 
by the six adjacent topological images of the massive particle, 
in a T 3 universe of three exactly comoving equal side lengths 
L a . The accelerations a3(r) are normalised to be constant if the 
dominant term is third order in r/L a [see Eq. (| 14f) ] and are 
shown against distance r < 0.01L a from the massive particle, 
where L = 3000/1" 1 Mpc. Positive accelerations are shown by 
plus symbols; negative accelerations are shown by circles. 

bution of accelerations. This distribution has a mean and 
standard deviation radial acceleration of 



(0.00 ±6.11 



Gm 
' i 2 " 



(13) 



and a skewness w 0.58(Gto/ L 2 )(r / L a ) 3 , where r = 
\/ x 1 + y 2 + z 2 . This expression is consistent with the lower 
scatterings of points in Figs. 7 and 8 of Roukema et al. 
(2007). Here, Fig. [T] shows a similar numerical calculation, 
using the radial component of the sum of 1/r 2 accelera- 
tions for the eight adjacent images, calculated directly from 
Eq. (10) of Roukema et al. (2007) without using the Taylor 
expansion. The acceleration is shown as 



a 3 (r) = 



Gm/Ll \L 



-3 



(14) 



for convenience. 

The analytical and numerical calculations are clearly 
consistent in showing that in the case of perfectly equal 
fundamental side lengths of a T 3 model, there is indeed 
a residual topological gravity effect. The residual accelera- 
tion is anisotropic, depending on the relation between the 
test particle's displacement from the massive particle and 
the orientation of the fundamental directions. The positive 
skewness implies that even though the mean radial accel- 
eration is close to zero, the mode and median are negative, 
i.e., test particles placed in random directions relative to 
the massive particle are more likely to be subject to a radial 
deceleration rather than a radial acceleration. The ampli- 
tude of the effect is third order, i.e., about (r/L a ) 2 times 
smaller than the residual effect that occurs for slightly un- 
equal fundamental lengths. If an upper estimate for r/L a 
at the present epoch is used, then this is a factor of about 
a million. 
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3.2. One pair of opposite topological images in S 3 /Y 

Before considering the full effect from a layer of topological 
images in the spherical cases, let us first consider the effect 
of just one pair of opposite topological images of the "local" 
massive particle. This is somewhat similar to the T x R 2 
case considered in Sect. 3.1.1 of Roukema et al. (2007), 
where the test particle lies along the geodesic joining the 
three images of the massive particle to one another in one of 
the spherical spaces, i.e., S 3 /Y for an appropriate holonomy 
group r. This may help us to understand the full sum from 
all the adjacent images. Similarly to Eq. (2) of Roukema 
et al. (2007), we can use scalar quantities. The acceleration 
to fifth order in t/Rq for the Poincare space is3 



G 



Rl 







-2 




















. sin (f 






sin (f - 


h i). 

















= G 



m f4cos(7r/5) r 



Rl 



sin 3 (7r/5) R C 
16 cos(7r/5) sin 2 (7r/5) 



3sin b (7r/5) 



24cos 3 (tt/5) / r V 



m 



15 - 9 fe) +sl - 6 (i)' 



3104 (£)' 



(15) 



This is similar to the flat case. The residual acceleration is 
proportional to the displacement to first order, which again 
is behaviour similar to that of a cosmological constant: 



r oc r. 



(16) 



The octahedral space has adjacent images at (tt/3)Rc, so 
the corresponding calculation for a single pair of images is 



G 



R 2 C 



Rc) 



16-\/3 „ m 

7 G —7T- 



G 



r 

R^ R~c 



2| — ) +— ( — ) 

Rc 5 \R C 



Re) 



(17) 



Equation (fTT)) is also valid for the adjacent images at the 
eight truncated corners of the fundamental domain in a 
truncated cube space, which are at (n/3)Rc- The other 
six adjacent images in the truncated cube space are at 
(n/4)R c , giving 



m 



Hi 



26 



7 (i 



65.1 



Rc, 



(18) 



Each of these calculations is for just one pair of topo- 
logical images adjacent to the "local" copy of the massive 
particle, lying in opposite directions. In the flat case, this 
can be thought of as a T x R 2 model, which is a valid 3- 
manifold. However, the spherical manifolds are not vector 
spaces, so a similar interpretation of the single-pair calcu- 
lation is not possible. 

On the other hand, Eqs ([To]). (fTTj) , and (JTHJ) hint at the 
form that numerical estimates for the sum of the weak-limit 
residual gravitational effect from all the adjacent images 
may take for a given spherical manifold. For S 3 /I*, the full 
set of adjacent images of the massive particle consists of 
six pairs of images. For a test particle displaced slightly 
in a random direction from the massive particle, the two 
images in a pair will be seen in nearly, although not ex- 
actly, opposite directions in the tangent 3-space at the test 
particle's location at r £ S 3 . The modulus of the (vector) 
residual acceleration induced by the nearly opposite pair 
should not be very different from the expression given in 
Eq. (fT5)) . although its expression using elementary algebra 
might appear complicated. 

This suggests that the scalar amplitude of the vector 
sum of all twelve accelerations is likely to involve terms 
of first, third, and fifth order in r/Rc, with coefficients 
of the order of magnitude of those in Eq. (|15p . However, 
this argument is not exact. A test particle displaced from 
the massive particle in an arbitrary direction does not, in 
general, lie along the great circle defined by a given pair 
of opposite images, and can at most lie along only one of 
the great circles defined by the six pairs of opposite images. 
Hence, terms with even powers of t/Rq could also appear 
in the Taylor expansion. 

These single-pair calculations might also be of interest 
for ill-proportioned (Weeks et al. 2004) positively curved 
spaces, e.g., the lens spaces L(p, q), with p,q G Z relatively 
prime, where p > 1 (e.g. Sect. 4, Gausmann et al. 2001). 
Since these spaces are not globally homogeneous, deriva- 
tions similar to those in Eqs (fl"5|) , (fTT|) . and (TT8"j) would 
be strictly valid only for points lying along the symmetry 
axis joining the centres of the two faces of the fundamen- 
tal dihedron (lens), i.e., where adjacent topological images 
are separated by a spatial geodesic of length 2irRc jp. This 
direction would therefore be unstable to a linear-term ac- 
celeration effect that would tend to expand it faster than 
other directions. This is qualitatively similar to the effect 
found in Sect. 3 of Roukema et al. (2007), according to 
which the residual acceleration would tend to equalise the 
three fundamental lengths of a T 3 model of slightly unequal 
fundamental lengths. 
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3.3. Well-proportioned spherical spaces: S 3 /T*, S 3 fO* 

s 3 /r 

As described in Sections 



and 



and 12. 3[ for a test particle at r 
in S 3 /T\ S 3 /0\ or S 3 /F, let us set R hm = l.l(ir/3)R c , 
Rhor = l.l(7r/3)i?c, or i?hor = 1. 1 (7r/5)i?c , respectively. 
Equation © is evaluated numerically for 200 randomly 
(logarithmically) distributed test particles at distances of 
up to 250/i -1 Mpc from the massive particle. Figures [H [21 
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Fig. 2. Residual acceleration (radial component) induced on a 
test particle near a massive particle by the eight adjacent topo- 
logical images of the massive particle, in a universe whose 3- 
manifold of comoving space is the octahedral space S 3 /T* . The 
accelerations as(r) are normalised to be constant if the domi- 
nant term is third order in r/L a [see Eq. (fT9)l ] and are shown 
against distance r < 250h~ 1 Mpc = 0.01-Rc from the massive 
particle. Positive accelerations are shown by plus symbols; neg- 
ative accelerations are shown by circles. 
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Fig. 3. As for Fig. [2] for the truncated cube space S 3 /O* . 

and [4] show the residual accelerations, scaled by (r/Rc) 3 
so that they should be approximately constant if the dom- 
inating term in r/Rc is the third order term, i.e., 



o 3 (r) 



Gm/R 2 c \R, 



is shown. 

It is clear that both the octahedral and truncated cube 
spaces have similar Taylor series behaviour to that of T 3 . 
The linear term cancels, but a third-order-dominated resid- 
ual acceleration remains. On the other hand, it is clear from 
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Fig. 4. Residual acceleration (radial component), as for Fig. [2] 
for the Poincare dodecahedral space S /I*. 




Fig. 5. Orthogonal residual acceleration for the Poincare do- 
decahedral space S 3 /I*, as for Fig. [3] shown using 50-bit (x), 
60-bit (o) and 70-bit (+) precision in the significands of the 
floating-point numbers. 



Fig. [4] that the Poincare dodecahedral space has a residual 
acceleration that is much weaker than those of the other 
three spaces, and that this is poorly modelled as a third or- 
der term in r/Rc . Instead, the constant slope of the relation 
in Fig. 2] strongly suggests that the residual acceleration for 
the Poincare space is dominated by a fifth order term, i.e., 



(19) 05 (r) 



Gm/R 2 c 



r 



-5 



(20) 



is approximately constant and terms lower than the fifth 
order cancel. 

The remaining orthogonal component of the residual 
acceleration is of a similar order of magnitude to that of 
the radial component. Figure shows that provided that 
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Table 1. Statistical characteristics of coefficients at of the dom- 
inant (i-th order) term in the radial and orthogonal compo- 
nents of the residual acceleration r in perfectly regular well- 
proportioned spaces, for approximately isotropic displacements 
r. a 



space 


term 1 




II /+° 




a U) 




t! 


3-torus 


(r/L a ) 


3 


II 


0.00 


0.00 


6.11 


0.58 


3-torus 


(r/La) 


3 


± 


6.39 


0.00 


6.82 


-0.42 


octahedral 


(r/Rc, 


) 3 


II 


0.01 


0.01 


3.18 


-0.58 


octahedral 


(r/Rc. 


f 


± 


3.33 


0.00 


1.26 


-0.42 


tr.cube 


(r/Rc. 


f 


II 


-0.05 


0.05 


14.24 


0.58 


tr.cube 


(r/Rc. 


I 3 


_L 


14.94 


0.02 


5.63 


-0.42 


dodec/num g 


(r/Rc, 


I s 


II 


-0.30 


0.53 


288.29 


0.75 


dodec/alg h 


(r/Rc 


) S 


II 


0.00 


0.01 


288.26 


0.74 


dodec/num g 


(r/Rc, 


) S 


_L 


286.33 


0.22 


121.65 


-0.37 



a coefficients a; as defined in Eqs (|14p , (|19p , and (|20|l ; these are 
approximately constant with respect to r; the constant factor of 
Gm/l? a for the T 3 model or Gm/R^ for the other models has 
been ignored here; all values shown are dimensionless 
b dominant i-th power of displacement, as derived in this paper 
c radial || or orthogonal _L component 

d standard error in the mean <T/ a > = ffi/y/N — 1 for iV 3> 1 test 
particles 

° sample standard deviation 

f sample skewness 7; = ([(«i — ( a i))/ <T i] 3 ) 

g from 70-bit significand numerical calculations using Eq. 

h using the algebraic expression in Eq. (|2ip 



the component is calculated with at least 70-bit precision 
in the significand, it is dominated by the (r/Rc) 5 term. 

The previous figures were calculated using 70-bit signif- 
icand multi-precision floating-point operations. The effect 
of reducing the precision to 50 or 60 bits is clear in Fig. [5l 
i.e., for the orthogonal component of the residual accelera- 
tion in the Poincare space. Noise created by the precision 
limit enters the calculation, for example, when converting 
from a double-precision position in R of the test particle 
to a multi-precision position. Apart from special cases, the 
finite precision representation of a typical, pseudo-random 
particle position in R 4 will place the particle at a position 
slightly offset from the physical 3-surface S 3 , because of the 
limited precision. If the particle is not located exactly on 
S 3 , then this induces an acceleration that is orthogonal to 
the tangent 3-plane at r. The numerical noise component 
appears on the left of Fig. [5] as a constant value of r, i.e., 
3 (r) oc (r/R c )~ 3 . 

3.3.1. Analytical estimate for the Poincare space 

Analytical evaluation of Eq. ([5]) to fifth order for 
the Poincare dodecahedral space using a computer 
algebra system__ confirms the numerical estimates 
shown in Figs. 2] and [5] The highest order residual 
acceleration from the adjacent topological images is 



5 The script is available online at 

http : //adjani . astro .umk.pl/GPL/dodec/PDS_residual 
Version 1.0 was used in this paper. 
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12^(297^ + 655) / r \ 5 

125^5- V5 \ R c) 
| [70 y 4 + (42 VE + 70) x 2 y 2 - (14 VE + 70) y 2 

+ (21 VE - 7) x 4 - 28 VE x 2 + 7 VE + 5] x, 
[70 z 4 + (42 VE + 70) y 2 z 2 - (14 VE + 70) z 2 
+ (21 VE - 7) y 4 - 28 VE y 2 + 7 VE + 5] y, 
[70 x 4 + (42 VE + 70) x 2 z 2 - (14 VE + 70) x 2 
+ (21 VE - 7) z 4 - 28 VE z 2 + 7 VE + 5] z, 

0}, (21) 

where the massive particle is at (0,0,0,1), 
the test particle displaced from it in an arbi- 
trary direction by a small amount is located at 
[sin(r/Rc)x, sin(r/Rc)y, sin(r/Rc)z, cos(r/Rc)], and 
x 2 + y 2 + z 2 = 1. The highest order term in the fourth 
component of this expression is a sixth order term in 
powers of r/Rc- The fifth order exact expression can 
be written with approximate decimal coefficients as 

r re (r/R c f [(4304.3x 4 + 7540. 5y 4 + 17657.2xV 
-6744.5a; 2 - 10912.8y 2 + 2224.7)a;, 
(4304.3y 4 + 7540. 5z 4 + 17657.2y 2 z 2 
-6744.5y 2 - 10912.8z 2 + 2224.7)y, 
(4304.3z 4 + 7540.5a; 4 + 17657.2z 2 a; 2 
-6744.5z 2 - 10912.8a; 2 + 2224.7)z, 
0]. (22) 

The radial and orthogonal components of the residual 
acceleration can now be calculated as follows. Firstly, a 
small numerical value r/Rc = 10~ 3 is chosen in order to 
evaluate the unit tangent vector at r pointing towards the 
massive particle. The sign of the latter is reversed in order 
to point away from the massive particle. Secondly, pseudo- 
random numerical values for x, y, z are substituted into the 
full vector residual acceleration [Eq. (|2"Tj) or (|2"2"jl ] and the 
scalar product with the unit tangent vector is taken in order 
to obtain the radial component of the residual acceleration. 
This is subtracted from the full residual to obtain the or- 
thogonal component. 

3.3.2. Statistical description 

For an isotropic distribution of the directions of displace- 
ment of the test particle, Table [1] lists characteristic statis- 
tics of the radial and orthogonal components of the residual 
acceleration for the four well-proportioned spaces, primar- 
ily calculated from direct numerical estimates using Eq. §§§ . 
By construction, the orthogonal component is necessarily 
non-negative. Since estimation of the parameters of the 
Poincare space residual acceleration is sensitive to numer- 
ical precision limits, Table [T] also lists parameters for the 
Poincare space radial component estimated using Eq. (|21[) 
and the unit tangent vector described above. Within the 
uncertainties indicated by the standard error in the mean 
c/ ai \, the two estimates agree. 
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4. Discussion 

The results shown above are quite sensitive to small errors 
in the positions of the adjacent topological images. This 
provides a useful check of the calculations, since a small 
arbitrary error is most likely to yield a stronger residual 
acceleration than a weaker one. Increasing or decreasing 
the distances from the massive particle to the two mem- 
bers of one of the opposite pairs of topological images by 
0.1% is sufficient to destabilise the nearly perfect equilib- 
rium defined by the full set of adjacent topological images 
for a given space. The residual acceleration reverts to being 
dominated by (r/Rc) 1 rather than (r/Rc) 3 or (r/Rc) 5 - 

This sensitivity to small changes in the distances to 
topological images is physically interesting. As found in 
Sect. 3 of Roukema et al. (2007), the residual acceleration in 
a T 3 model of slightly unequal fundamental lengths is domi- 
nated by a linear term in r/L a , which tends to equalise the 
three fundamental lengths. The same relation applies for 
the spherical models. These are very well balanced in the 
sense that the linear term cancels, and in the case of the 
Poincare space, the cubic term cancels too, provided that 
these spaces are perfectly homogeneous. A small decrease 
in the injectivity diameter in one direction implies stronger 
acceleration for a test particle, favouring a return towards 
perfect isotropy in the fundamental lengths. Conversely, a 
small increase in the injectivity diameter implies weaker ac- 
celeration, again favouring a return towards isotropy. The 
residual acceleration tends to encourage the space to return 
to perfect "residual acceleration equilibrium" , in which the 
residual acceleration disappears down to the third or fifth 
order. 

Clearly, these effects are negligible in the present epoch 
at the level of observational detectability for the next few 
decades. Moreover, given the present empirical interest in 
the Poincare space, the prediction of the perfectly isotropic 
model would be an absence of a residual acceleration effect 
to a much higher observational accuracy than that for the 
effects predicted by the other three well-proportioned mod- 
els, which in turn would be more difficult to detect than the 
residual acceleration effect from ill-proportioned models. 

However, the Universe is certainly inhomogeneous. If 
evidence continues to accumulate for the Poincare space 
model with successively more accurate estimates of the 
model's astronomical coordinates, following the estimates 
made in Roukema et al. (2008a, 2008b), then several ob- 
servational approaches could be used to test the predic- 
tions of the residual acceleration effect. With a sufficient 
level of precision, matched circles and/or annuli on the sky 
should yield slightly different fundamental lengths along 
the six axes. The residual acceleration effect should consist 
of a positive or negative extra acceleration for shorter or 
longer fundamental lengths, respectively. Surveys of tracers 
of large-scale structure lying along the different geodesies 
to topological images of, for example, the Virgo cluster, 
could also be used to estimate which directions should have 
slightly longer than average geodesies and which should 
have slightly shorter than average geodesies. All three es- 
timates should agree with one another, provided that no 
other sources of random and systematic error interfere. 
Hence, the residual acceleration effect potentially offers a 
physical mechanism for testing a multiply-connected model 
of comoving space, as an alternative to relying on geometri- 



Table 2. Components of the residual acceleration r, as for 
Tablefl] for the truncated cube space, separated by fundamental 
domain face shape. 
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II /-L 


(Oi) 


a (o.i) 


o, 


Ji 


octagonal a 


(r/Rc) 3 


II 


-0.00 


0.04 


17.46 


0.58 


octagonal 


(r/Rc) 3 


_L 


18.27 


0.02 


6.89 


-0.43 


triangular 10 


(r/Rc) 3 


II 


0.00 


0.01 


3.19 


-0.58 


triangular 


(r/Rc) 3 


_L 


3.33 


0.00 


1.26 


-0.43 



a contribution from the octagonal (truncated square) faces 
contribution from the triangular faces 

cal effects (individual or statistical identification of multiply 
imaged objects or regions of space). 

Another interesting question is whether or not the bet- 
ter "balancing" of the Poincare space could have had any 
role in the evolution of a preferred topology in the early 
Universe, especially during the quantum epoch. Does the 
Poincare space occupy a dynamically more stable state 
than the other well-proportioned spaces, which in turn are 
more stable than ill-proportioned spaces? This would most 
likely require significant inhomogeneities in the very early 
Universe, which themselves might have been subsequently 
removed by the residual gravity acceleration effect itself. 

Apart from the possible physical consequences of this 
result, it may be interesting to ask if these results are rea- 
sonable from an intuitive mathematical point of view. It is 
useful to consider this in terms of pairs of opposite images, 
since if the test particle lies along the geodesic joining the 
local copy and the two opposite distant copies of the mas- 
sive particle, then the single pair yields a linear-dominated 
residual acceleration. In order for the linear (or higher or- 
der) terms in the residual acceleration to cancel, the accel- 
eration components that are approximately orthogonal to 
the geodesic joining the test particle to the two members of 
a given opposite pair need to be able to cancel the approxi- 
mately radial components of residual accelerations in other 
pairs of images. In the T 3 case, three nearly orthogonal 
pairs are sufficient for the linear terms to completely cancel 
in this way (provided that the fundamental domain is per- 
fectly regular). Given that S 3 /T*, S 3 /0* 1 and S 3 /I* have 
four, seven, and six opposite pairs, respectively, it is rea- 
sonable to imagine that the radial components of the linear 
residual-acceleration term of a given pair are cancelled by 
the orthogonal components of the other pairs, since there 
are many of them, distributed symmetrically. 

It is also possible to imagine that with a high enough 
number of opposite pairs and sufficient symmetry, a higher 
order term such as the third order might also disappear, 
by means of fine balancing between the different directions. 
The truncated cube space S 3 /0* and the Poincare space 
S 3 / 1* have the highest numbers of opposite pairs, or equiv- 
alently, the highest numbers of pairs of faces of their fun- 
damental domains. However, the truncated cube cannot be 
a regular polyhedron; its faces are six truncated squares 
and eight triangles. All the faces of the fundamental do- 
main of the Poincare dodecahedral space have the same 
shape: a pentagon. In this sense, it is more symmetrical 
than the truncated cube space. So requiring both a fun- 
damental domain with a high number of faces and a high 
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level of symmetry would seem to provide a qualitative pair 
of conditions to explain why the Poincare space is more 
finely balanced than the other well-proportioned spaces. 

However, could either the subset of topological images 
corresponding to only the octagonal faces of the fundamen- 
tal domain of the truncated cube, or the subset correspond- 
ing to only the triangular faces, nevertheless cancel to fifth 
order, and become hidden because the other subset rein- 
troduces a third order term? A calculation similar to those 
developed above shows that both subsets separately retain 
the third order as the dominant order. Table gives sta- 
tistical characteristics of the third order coefficients for the 
two subsets. So the Poincare space retains its uniqueness in 
cancelling down to fifth order. 

An aspect of the residual acceleration as expressed in 
Eqs (|2"Tj) and (f2"2"| that may initially seem counterintuitive 
is the permutational symmetry between the x, y, and z com- 
ponents, since a dodecahedron is not always thought of as 
having cubical symmetry. However, an appropriate orien- 
tation of the dodecahedron in R 3 , such as that used in our 
computer algebra script, shows an (x, y, z) symmetry. The 
twelve adjacent topological images consist of three quadru- 
plets orthogonally projected to R 3 , [0, ±1/2, ±( v / 5 - l)/4], 
[±(V5-l)/4,0,±l/2], and [±1/2, ±{y/E- l)/4, 0]. Each of 
these quadruplets lies in a 2-plane in R 3 , i.e., x — 0, y = 0, 
and z = 0, respectively. In this orientation, the "top" and 
"bottom" of the dodecahedron are edges, not faces. These 
two edges, along with four others, can be used to inscribe 
the dodecahedron in a cube (in flat space). 

5. Conclusions 

In Roukema et al. (2007), it was found that a residual 
gravity acceleration effect exists for T x R 2 but cancels 
to first order for a perfectly regular T 3 model. Here, it 
was found that the residual acceleration in T 3 does not 
completely cancel. In T 3 and two other well-proportioned 
spaces, the octahedral space S 3 /T* and the truncated cube 
space S 3 fO*, the residual acceleration exists as a third or- 
der effect in the Taylor expansions in r/L a and r/Rc re- 
spectively. A reasonable upper limit to r/L a = r/Rc can 
be set by considering test particles in voids of large-scale 
structure, i.e., at typically a few megaparsecs from the most 
massive nearby object, and an injectivity diameter for the 
3-manifold of comoving space of at least a few gigaparsecs, 
i.e., r/L a = r/Rc ~ 10~ 3 . Hence, at the present epoch, 
these three well-proportioned spaces are about a million 
times better balanced by this dynamical criterion than ill- 
proportioned spaces. 

The Poincare dodecahedral space S 3 /I*, presently a 
candidate 3-manifold for comoving space favoured by sev- 
eral observational analyses, has been found to be even more 
exceptional. Its residual acceleration is dominated by the 
fifth order term of amplitude - ±300(Gm/ R c )(r / R c ) 5 
[Figs. SI TableEQ Eqs jH]), ([22])]. This makes it about ten 
thousand times better balanced than the other three well- 
proportioned spaces, i.e., about 10 10 times better balanced 
than ill-proportioned spaces. Moreover, perturbations to 
this equilibrium favour a return to the equilibrium state. 
Are these clues towards a theory of cosmic topology? 
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Appendix A: Residual gravity in the Poincare 
space: computer algebra script 

The following script evaluates the residual acceleration on 
a test particle near a massive particle due to the adja- 
cent topological images of the massive particle, as given 
in Eq. ©, for the Poincare dodecahedral space. The script 
is written using the syntax of the computer algebra system 
MAXIMA, version 5.10.0. 

/* 

# topograv - numerical tests of residual gravity 

# force induced by spatial topology 
# 

# Copyright (C) 2009 Boud Roukema 
# 

# This program is free software; you can 
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# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 
# 

*/ 



redistribute it and/or modify it under the 
terms of the GNU General Public License as 
published by the Free Software Foundation; 
either version 2, or (at your option) any 
later version. 

This program is distributed in the hope that 
it will be useful, but WITHOUT ANY WARRANTY; 
without even the implied warranty of 
MERCHANTABILITY or FITNESS FDR A PARTICULAR 
PURPOSE. See the GNU General Public License 
for more details. 

You should have received a copy of the GNU 
General Public License along with this 
program; if not, see: 

http : //www . gnu . org/licenses/gpl . html 



/* avoid imaginary sqrt ambiguities */ 
radexpand: false; 

normalise_v4(x) := [ x [1] /sqrt (x . x) , 
x[2] /sqrt(x.x) , x [3] /sqrt (x . x) , 
x[4]/sqrt(x.x) ] ; 



/* normalised tangent at 4-vector a pointing 

towards 4-vec b */ 
tangent_S3_at_a(a, b) := 
normalise_v4(b - ((a . b) * a) ) ; 

/* spherical-newtonian acceleration at a from 

image at b */ 
accel_at_a(a,b) := tangent_S3_at_a(a, b) / 
(1- (a . b)~2); 

/* zeroth image of massive particle is at 
[0,0,0,1] */ 

/* number related to the Golden Ratio */ 
phi_ml4 : (sqrt (5) -1) /4 ; 

phi_pl4 : (sqrt(5) + l)/4 ; /* = cos(°/„pi/5) */ 

/* 12 adjacent images of massive particle; 
modulus must be 1 since the physical space 
is S~3, not R~4 */ 



image 1 
image2 
image3 
image4 
image5 
image6 

image7 
image8 
image9 
image 10 
image 11 
image 12 



[ , 1/2, -phi_ml4, phi_pl4 ] 

[ , -1/2, -phi_ml4, phi_pl4 ] 

[ -1/2, phi_ml4, 0, phi_pl4 ] 

[ 1/2, phi_ml4, 0, phi_pl4 ] 

[ -phi_ml4, 0, -1/2, phi_pl4 ] 

[ -phi_ml4, 0, 1/2, phi_pl4 ] 

[ , -1/2, phi_ml4, phi_pl4 ] 

[ , 1/2, phi_ml4, phi_pl4 ] 

[ 1/2, -phi_ml4, 0, phi_pl4 ] 

[ -1/2, -phi_ml4, 0, phi_pl4 ] 

[ phi_ml4, 0, 1/2, phi_pl4 ] 

[ phi_ml4, 0, -1/2, phi_pl4 ] 



accel_f ull 
accel 
accel 
accel 
accel 
accel 
accel 
accel 
accel 
accel 
accel 
accel 
accel 



(vec_m) := 
_at_a(vec. 
_at_a(vec. 
_at_a(vec. 
_at_a(vec. 
_at_a(vec. 

at_a(vec. 

at_a(vec. 

at_a(vec. 
_at_a(vec. 
_at_a(vec. 
_at_a(vec. 

at_a(vec. 



m, image 1) 
m, image2) 
m, image3) 
m, image4) 
m, image5) 
m, image6) 
m, image7) 
m, image8) 
m, image9) 
m.imagelO) + 
m,imagell) + 
m, imagel2) ; 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



/* test particle near to zeroth image of 
massive particle, in arbitrary direction 
(x,y,z); modulus must be 1 since the physical 
space is S~3, not R~4 */ 
vec_m : [ sin(rr)*x, sin(rr)*y , sin(rr)*z, 
cos(rr) ] ; 



/* Taylor expand to fifth order */ 

tay : taylor (accel_f ull(vec_m) ,rr,0,5) ; 

/* simplify x, y components */ 

tay_x : gf actor (rectf orm(expand(tay [1] ))) , 

z : sqrt(l-x~2-y~2) ; 
tay_y : gf actor (rectf orm(expand(tay [2] ))) , 

x : sqrt(l-y~2-z~2) ; 
/* simplify z component; to retain symmetry, 
reverse alphabetical order of (x,z) 
by temporarily substituting (b,a); 

*/ 

subst(a,z,subst(b,x,tay[3] )) ; 
gf actor (rectf orm(expand(°/ ) ) ) , 

y : sqrt(l-a~2-b~2) ; 
tay_z : subst(z,a,subst(x,b,°/ )) ; 
/* simplify w component */ 
tay_w : expand (tay [4] ) , 

x : sqrt(l-y~2-z~2) ; 

/* still fully algebraic */ 
residual_fully_alg : 

[tay_x, tay_y, tay_z, tay_w ]; 

/* 212-bit precision for big floats lbO etc. */ 
fpprec : 64; ratepsilon : lb-64; 

/* numerical coefficients in front of x,y,z */ 
residual_xyz_alg : residual_fully_alg, bfloat; 

/* put some weakly pseudo-random numbers into 
x,y,z, while requiring the R~3 modulus to 
be 1 so that vec_m is on S~3; */ 

x: 2b0*random(l .0)-lb0; 

theta: 2b0*°/„pi *random(l . 0) ; 

y: sin(acos(x)) * cos(theta), bfloat; 

z: sin(acos(x)) * sin(theta) , bfloat; 

residual : residual_xyz_alg, bfloat; 



/* sum of accelerations */ 



/* components */ 
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radial_n: -tangent_S3_at_a(vec_m, [0, 0, 0, 1] ) , 
rr : lb-3; 

residual_radial : expand (residual . radial_n) ; 

ortho : expand (residual - 

residual_radial * radial_n) ; 
residual_ortho : expand (sqrt (ortho . ortho) ) ; 

/* result summary */ 

print ( "vector residual (coeff, x,y,z algeb.)"); 
residual_f ully_alg ; 

print ( "vector residual (x,y,z algebraic)"); 
residual_xyz_alg; 

print ( "vector residual (x,y,z numeric)"); 
residual ; 

print ( "radial residual (x,y,z numeric)"); 
residual_radial ; 

print ( "orthogonal residual (x,y,z numeric)"); 
residual_ortho; 



